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ABSTRACT 

Forecasting large Solar Energetic Particle (SEP) events associated with shocks driven by fast coronal 
mass ejections (CME) pose a major difficulty in the field of Space Weather. Besides issues associated 
with CME initiation, the SEP intensities are difficult to predict, spanning 3 orders of magnitude at 
any given CME speed. Many lines of indirect evidence point to the pre-existence of suprathermal 
seed particles for injection into the acceleration process as a key ingredient limiting the SEP intensity 
of a given event. This paper outlines the observational and theoretical basis for the inference that 
a suprathermal particle population is present prior to large SEP events, explores various scenarios 
for generating seed particles and their observational signatures, and explains how such suprathermals 
could be detected through measuring the wings of the H I Ly-a line. 

Subject headings: Acceleration of Particles - Shock Waves - Sun: corona - Sun: coronal mass ejections 
(CMEs) 



1. OVERVIEW 

An extensive body of evidence identifies shocks driven 
by very fast coronal mass ejections (CMEs) within a few 
solar radii of the Sun as the primary particle accelerators 
in large, gradual solar energetic particle events. The vari- 
ability of SEP events and the wide range of particle and 
contextual information offer a powerful tool for exploring 
shock acceleration, which is ubiquitous in astrophysical 
plasmas. A predictive capability for SEPs would also be 
of great applied interest. Diverse lines of evidence indi- 
cate that rapid production of large intensities of high- 
energy particles is greatly enhanced when the pre-event 
environment has been primed with a population of ions 
having energies well above the typical thermal particle 
energy, usually in the range from a few to tens of keV in 
the solar corona. However, at present we have no direct 
evidence that such "suprathermal" ions actually exist in 
the corona in numbers sufficient to serve as seed parti- 
cles for diffusive shock acceleration (DSA). Instead, all 
evidence for suprathermals in the corona is inferential. 

A primary science objective must be to detect 
suprathermal protons in the corona at altitudes where 
the production of high-energy SEPs begins. We discuss 
the use of high-throughput UV spectroscopy to quan- 
tify broadening of the Ly-a line resulting from charge 
exchange between suprathermal protons, with energies 
in the range 1-5 keV, and neutral hydrogen atoms in 
the corona. This method provides remote-sensing of the 
coronal conditions that allow large SEP events to oc- 
cur. Concentrating our efforts on the seed particle prob- 
lem, we envisage a robust, simplified design optimized 
to detect coronal suprathermals, with an effective area 
at H I Ly-a being more than 500 x larger than that 
of instruments previously flown. Such instrumentation 
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would also provide a prototype for future space-weather 
instrumentation, and offers great synergy with the He- 
liophysics System Observatory. By determining the den- 
sity of suprathermal protons in the solar corona and the 
variations thereof, the connection between suprathermals 
and SEP fiux, and the physical processes that create 
suprathermals could be investigated directly. 

In this paper we first outline the current indirect evi- 
dence for suprathermal seed particles, both observational 
and theoretical, in section 2. In section 3, we examine re- 
quirements on the seed particle distribution function. We 
introduce the "kappa" distribution to model wave gener- 
ation ahead of a CME shock, and radiation transfer ef- 
fects on the Lyman a line profile by suprathermal atoms 
in the corona. The kappa distribution need not be the 
exact description of the suprathermal tail, but it provides 
an easily calculable framework while capturing important 
features of the likely seed particle distribution. Section 4 
explores some of the possible mechanisms for generating 
seed particles, the different observational signatures they 
would have, and the extent to which they can meet the 
requirements already discussed in section 3. In section 
5, we describe some of the secondary science, concerned 
principally with electron acceleration, that arises. Sec- 
tion 6 concludes with considerations of the importance of 
such science to space weather and the prediction of SEP 
events, and the experimental implementation of such a 
strategy. 

2. INTRODUCTION: THE SEP EFFECTIVENESS OF CMES 

At any given CME-speed, the observed SEP intensi- 
ties span more than three orders of magnitude (Reames 
2000). There are many strong observational reasons to 
think that an important factor in this variability other 
than CME speed, is variation in the available suprather- 
mal seed population. The intensity of remnant energetic 
particles from preceding events is the only factor discov- 
ered so far that helps organize event-to-event variabil- 
ity in the peak SEP intensities at 1 AU (Kahler et al. 
1999; Kahler 2001). Gopalswamy et al. (2004) found 
that large SEP events were generally preceded by an- 



2 



other CME erupting from the same active region within 
the preceding 24 hours (see also Kahler & Vourhdas 
2005). Chver (2006) found that an enhanced pre-event 
background also favored the occurance of Ground Level 
Enhancements (GLEs). The previous CME had appar- 
ently "pre-conditioned" the corona to foster the produc- 
tion of a large SEP event. Mewaldt et al. (2012) found 
that an upper-limit on the fluence of SEP Fe at 10- 40 
MeV/nucleon can be predicted from the Fe intensity ob- 
served at Earth at 0.01-2.0 MeV/nucleon on the previ- 
ous day. Mason et al. (2005) reported a 100-day sur- 
vey at 1 AU of variation in the intensity of suprather- 
mal Fe ions at ~30 keV/nucleon, compared to a simul- 
taneous survey of the bulk solar-wind Fe. Whereas the 
bulk solar- wind density varies by only a factor of ten, the 
suprathermal intensities span three orders of magnitude. 
Thus, the suprathermal variability is of the same magni- 
tude as the scatter in SEP intensity variations as seen 
by Reames (2000). But, of course, these suprathermals 
were observed at 1 AU, not in the corona. A number 
of observational studies have also identified suprather- 
mals from flares, which reveal themselves by distinctive 
compositional signatures (i.e. high Fe/0 and •^He/'*He 
compared to solar- wind and coronal values), as another 
significant component of the shocks's seed population in 
large gradual SEP events (Mason et al. 1999; Tylka et 
al. 2001: 2005; Desai et al. 2006ab; Tylka & Lee 2006a; 
Sandroos & Vainio 2007). In fact, evolving shock geome- 
try, coupled with a compound seed population compris- 
ing suprathermals from the corona, previous CMEs, and 
from flares, can account for many otherwise unexplained 
aspects of compositional and spectral variability in large 
SEP events (Tylka & Lee 2006ab; Sandroos & Vainio 
2007); 

Tvl ka fc Led (|2006b[ ) demonstrate that the injection 
profile into the heliosphere at the s un f or ^ GeV protons 
, as deduced bv lBieber et"al] (|2004[ ) and lSaiz et all (|2005[ ) 
from time-intensity profiles and angular anisotropies in 
ground level events from the "Spaceship Earth" world- 
wide network of neutron monitors, occurs when the CME 
shock is within a few (2-3) solar radii. Other studies, 
based on satellite observations of velocity dispersion in 
the first arrival of particles at earth, yield times of par- 
ticle departure from the sun that agree wi th the neu- 
tron monitor results to within a few minutes () Tvlka et al.l 
120031 : iReamesI l2009| ) . If we are to discover how shocks 
produce these very high energy particles, this is the re- 
gion of the corona where we must understand not only 
the shocks but also the characteristics and contents of 
the medium through which they move. 

On the theoretical side, there are also good reasons to 
believe that a pre-existing suprathermal seed population 
can account for the variability of CME SEP characteris- 
tics. One of the long-standing problems in the theory of 
DSA is how to inject particles from the thermal popula- 
tion into the acceleration process (e.g. Malkov & Drury 
2001, Malkov & Voelk 1995). Because DSA involves par- 
ticles being able to scatter across the shock there is a 
minimum energy below which a particle cannot be ac- 
celerated by DSA. According to Zank et al. (2006), for 
fast CME driven shocks within 1 A.U. (as modeled in 
Zank et al. 2000), this is of order lO^eV for parallel 
shocks and 20 times higher for perpendicular shocks, but 



can of course vary with shock parameters. Many authors 
have explored how one can achieve DSA through "ther- 
mal leakag e" from the thermal popu l ation (e.g. Ellison et 
al. 1995). IGahnskv fc Shevchenkd (poTlh give a model 
for proton injection directly from the upstream medium 
at oblique shocks, and apply it to solar wind events. We 
will discuss this in more detail below, but remark for 
now that such ideas may have more currency in connec- 
tion to astrophysical shocks. However, if a population of 
suprathermal particles already exists, as suggested by ob- 
servation of elemental abundances in gradual SEP events 
(references cited above), then the need to energized from 
the solar-wind thermal population (the injection prob- 
lem) becomes less relevant or non-existent. 

This is demonstrated in Fig. 1, by new calculations 
based on work in Ng & Reames (2008). Proton accel- 
eration at a parallel shock is simulated over the first 10 
minutes of CME motion outward from 3 i?©, for two dif- 
ferent seed particle populations (red and blue lines on 
the left hand panel). The middle and right hand panels 
show the ensuing proton spectrum at different times for 
the lower and greater seed particle density cases respec- 
tively. The nonlinear aspect of shock acceleration means 
that an order of magnitude difference in seed particle 
density translates into many orders of magnitude dif- 
ference in SEP densities at energies above the Coulomb 
barrier; 10 - 20 MeV. Consequently, very different radi- 
ation hazards would be expected in these two cases. For 
quasi-perpendicular shocks, where the injection thresh- 
old for efficient acceleration at the shock is higher (Zank 
et al. 2006), the need for suprathermal seed particles 
and the differing resulting radiation signatures would be 
e ven more a. c ute. 

IGiacalond (|2005allbh argue that a seed particle dis- 
tribution is unnecessary, even for perpendicular shocks, 
in the case tha t the shock is run ning into an irregular 
magnetic field. iGiacaloni ()2005aD considers high Mach 
number shocks running into magnetic field with fiuc- 
tuations 6B/ B between 0.1 and 1 using a test parti- 
cle approach. IGiacalond ()2005br ) studies a lower Mach 
number shock {Ma = 4, closer to that appropriate for 
CMEs under consideration here) using hybrid simula- 
tions, with similar results. It might be argued that 
these works are sim ply substituting "seed turbulence" 
for "seed particles". iGargate &: Spitkovskvl (j2012i) pro- 
vide a series of hybrid simulations of shocks with differing 
obliquity propagating into initially undisturbed cold up- 
stream media. They find particle acceleration efficiency 
decreasing with increasing shock obliquity, with essen- 
tially no acceleration occurring at strictly perpendicu- 
lar shocks. Observations of suprathermal protons (47 - 
65 keV) associated with shocks in the solar wind at 1 
A.U. show no si gnificant dependence on shock obliquity 
(jCiacalond 12012) . at least for the Ma > 3 shocks in- 
cluded in the sample. Analysis of suprathermal protons 
observed at a set of quasi-parallel solar wind shocks, also 
at 1 A.U. (jNeergaard Parker &: Zankll2012l ) revealed in- 
jection energies in the range 1-3 keV. These authors con- 
cluded that injection proceeds directly from the upstream 
Maxwellian distribution, and that a separate population 
of "seed" particles is not necessary. However we empha- 
size that both these works concern shocks observed at 
1 A.U., whereas in this paper we focus on shocks much 
closer to the Sun. 



3 



Although theoretical opinion is far from unanimous, it 
is generally recognized that the presence of a suprather- 
mal population greatly facilitates the effects of shock 
acceleration. Both observations and theoretical studies 
suggest that the availability of coronal suprathermal ions 
is likely to be a major factor in SEP event-to-event vari- 
ability. But this notion must necessarily be regarded 
as merely an attractive hypothesis until we find direct 
evidence for coronal suprathermals, independent of the 
SEPs themselves. 

3. OBSERVATIONAL SIGNATURES OF SEED PARTICLES 

3.1. The "Kappa" Distribution 

Many authors (e.g. iZank et"al]|2006[ ) discuss injection 
requirements in terms of a minimum particle energy. 
Spectroscopically, what we actually observe, for instance 
in a line profile, is a velocity distribution function of ions, 
and not an energy of a specifically suprathermal popula- 
tion. Here we discuss the description of a suprathermal 
ion population by a "kappa" distribution. Normalized 
over three dimensional velocity space, the function is il- 
lustrated in Figure 2 and is given by 
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where F (k) is the Gamma function of argument k. As 
K — > cx), the distribution tends to a Maxwellian with 
thermal speed vth- For lower values of k, the distribu- 
tion develops stronger wings; i.e. has more particles at 
higher velocities than the corresponding Maxwellian. If 
the integration over velocity extends to infinity, we re- 
quire K > 3/2 to keep the number density finite, and 
K > 5/2 for finite energy density. Lower values of k 
require a high-velocity cutoff to keep the suprathermal 
particle number and energy densities finite. 

In the same way that the Maxwellian distribution can 
be derived from the Boltzmann equation in terms of 
the coefficients of dynamical friction and velocity diffu- 
sion, the kappa distribution is derived with the inclu- 
sion of a v elocity diffusion coeffic ient associated with 
turbulence ([Hasegawa et al.l I1985D . This is done for 
electrons in lower-hybr id turbulence in Appendix B of 
iLaming Sz Lepril (I20p7fl. and can similarly be done for 
suprathermal ions. iTreumann &: Jaroschekl (12008) give 
more universal derivations based on an adaptation of the 
Gibb's distribution (a generalized Lorentzian) appropri- 
ate for long range correlations between particles in coUi- 
sionless plasmas. Equation 1 is similar to the generalized 
Loren tzian distribution given by equation 15 of lCranmeri 
(|1998f l. and is exactly the same if we identify k = kc + 1 
and 2kv\,j^ = {nc) wj^^, where kc and w^c are the value 
of kappa and the thermal width in Cranmer's expression 
(we have added the subscript C). Thus we are able to 
employ the redistribution functions for radiation transfer 
given in this reference.^ 

Many authors (e.g. iZank et"aLl |2006[) argue that in- 
jection energies for seed particles at quasi-perpendicular 
shocks should be significantly higher than is the case at 
quasi-parallel shocks. This suggests that ion distribu- 
tions deviating more strongly from Maxwellian distribu- 
tions (i.e. with a lower value of k giving more pronounced 
"wings" and a higher fraction of particles at higher en- 
ergies) should be present at quasi-perpendicular shocks. 



However it is difficult to estimate what the distribution, 
or more quantitatively the value of k, should be from such 
calculations. Here we adopt a criterion that the seed par- 
ticle distribution must be capable of generating parallel 
propagating Alfven waves when streaming upstream of 
a shock. This represents a minimum criterion for diffu- 
sive shock acceleration, in that upstream particles must 
be able to self-generate waves, so they can be scattered 
back downstream. Downstream of course, the shock it- 
self can generate turbulence, even though the particles 
may also do so. We envisage a population of seed parti- 
cles that is initially stationary in the upstream medium, 
being subsequently reflected further upstream upon in- 
teraction with a shock. Writing the angle between the 
shock velocity vector and the pre-shock magnetic field 
direction as 0;,„, we evaluate the u pstream grow th rate 
for parallel propagating waves fe.g. lMelroselll986l ) 



7=-///^-^Bin^a5(..-f7-fc||.||) 
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where a is the angle between the particle velocity vec- 
tor and the magnetic field, and other symbols have their 
usual meanings; u> is the wave frequency, J7 is the par- 
ticle cyclotron frequency, va is the Alfven speed, v and 
p are the particle velocity and momentum respectively, 
with subscripts _L and || indicating components perpen- 
dicular or parallel to the ambient magnetic field. We 
restrict our attention to parallel propagating wav es, as 
seen in the simulations of Gargatc & Spitkovskv t2012f ) 
and the in situ observations of iBamert et al.i ()2004i ) . In a 
coordinate system where the magnetic field defines the z 
axis, the x axis lies in the plane defined by the magnetic 
field and shock velocity vector, and the y axis is normal 
to this plane, the shock velocity vector has components 
{—Vs sin0f,n, 0, Vs cos^bn). The streaming solar energetic 
particle (SEP) distribution is then 



/ 



/o(l 



3va.v^ 



— /o (1 -I- 3^ (cos a COS 0b„ — sin a sin 9bn cos cf))) . (3) 

Substituting this into equation 2 and integrating over <j)^ 
we find 
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where is the number density of the streaming seed par- 
ticle distribution. Below, n will represent the stationary 
background ion density. For nonrelativistic cosmic rays, 
p = mv, V = r2/(fc||| cosa|) = VA^/{t^\cosa\) and so 
integrating over p 



4 



sin^ a. 



3— wcfl, cosa |cosa 



cos 6'6„ 



2kuj-^v^iJ cos^ a\ 



v^fl"^ sin^ a 3vcrv\Q sin^ a cos a 
v^f^uj'^ I cos^ a I 



cos dbn 



1 



2„2 



T/i 



cos^ a 



-(k.+ I) 



d(cos a). 



(5) 



Since the region of integration is — 1 < cos a < 1, the 
terms in cos a integrate to zero. Evaluation of the re- 
maining terms, either by analytic approximation or nu- 
merical integration, is aided by making the substitution 
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with -q = y/2KijjVTh/^VA, and using fi^/ojpj = and 
P = 2v^i^/v\, the integral becomes 
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We evaluate equation 7 numerically, assuming Q,/lo — 
max {M A cos 6 bm 5)- This comes from writing — 
k^^Vs cos9hn — ^MACOsOhn, and imposing a maximum 
value of w as 6* — )■ 7r/2, which is discussed further below. 
We take the limits of integration to be 7r/2 — arctanyy < 
X < 7^/2 + arctan?]. Analytically, 
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where rj' = sin tan ^ rj < 1. In the limit r; << 1 

- n 2/32 0.4 r (« _ 3/2) „^ COSWbn 

(9) 

suggesting k < 1.5Af^cos06„ (to keep 7 > 0) as a simple 
criterion for particle injection into the shock acceleration 
process at oblique shocks. A stronger non-thermal tail 
is required at higher shock obliquity, and at lower shock 
Alfven Mach number. The forgoing treatment is cor- 
rect for the case of a kappa distribution of ions flowing 
through a background Maxwellian, whose particle den- 
sity appears in the denominator of va- One can easily 
verify that all growth and damping rates tend to zero as 
K — 00 so the addition of a Maxwellian has no effect here. 
A streaming kappa distribution of density moving 
through a background of density n with the same kappa 
distribution gives k < 1.5Ma cosOhnn-r/ (nr + n), a more 



stringent requirement than the stationary Maxwellian re- 
quires. 

We have tacitly assumed << n, which need not 
be the case in reality. In this case the definition of 
VA, both in terms of the density entering the denomina- 
tor, and the frame of reference (stationary background 
or streaming ions) becomes ambiguous, and requires re- 
visiting the dispersion relations for MHD waves in the 
composite plasma. 

3.2. A Case Study: The 2005 January 20 Event 

To illustrate how wave growth ahead of CME shocks 
might evolve, and how the particle distributions respon- 
sible for such wave growth might be detected, we con- 
sider the 2005 January 20 event. The flare and asso- 
ciated CME and SEPs have been studied by many au- 
thors. iGrechnev et al.l (|2008[ ) give a survey of the flare 
evolution from imaging and timing data acquired by the 
Reuven Ramaty High Energy Solar Spectroscopic Imager 
(RHESSI; Lin et al. 2002) and the Transition Region 
and Coronal Explorer (TRACE; Handy et al. 1999), and 
study the propagating wave that develops into a shock 
using the Extreme-Ultraviolet Imaging Telescope (EIT; 
Delaboudiniere et al. 1995) 195 A channel and white 
light images from the Large- Angle Spectrometric Coron- 
(7^graph (LASCO; Brueckner et al. 1995) C2 and C3. El- 
em ent charge s t ates w ithin the CME ejecta are studied 
bv iLepri et al.l (|2012[ ). while othe r in situ ob servations 
(ia;are compreh ensively reviewed by iFouUon et a l. (2007). 
According to IGrechnev et al.l (|2008t ). the shock reached 
a velocity of about 2000 km s~^ at a heliocentric dis- 
tance of 5 i?0, and quite possibly ac hieved this speed at 
lower, less well observed, altitudes. iLepri et al.l (|2012j ) 
model the behavior of the Alfven speed with altitude. 
The CME wave first emerges as a shock at a heliocentric 
radius of about 1.3 i?©, where the Alfven speed is around 
2000 km s~^. The Alfven speed declines monatonically 
as the shock moves out, reaching Alfven Mach numbers 
o f 2 at about l.S Rp, an d 3 at 2.1 Rq. 

IGrechnev et"al] (|20 0g ) argue that all the electromag- 
netic emissions they analyze come from the flare site and, 
interestingly, also suggest that the flare was the location 
of the associated particle acceleration. We will take the 
view advocated above that SEPs are accelerated at the 
shock, aided by the existence of seed particles in the up- 
stream medium, and calculate the distribution function 
(i.e. the value of n) these seed particles must have in 
order to initiate particle acceleration. The 2005 Jan- 
uary 20 event was the last in a series of four CMEs from 
the same active region occurring on January 15, 17, 19, 
and 20, and had by far the largest SEP productivity. 
IGrechnev et al. ' (2008) comment that "seed populations 
from previous CMEs do not appear to be drastically dif- 
ferent for the 17-20 January events", but in the absence 
of remotely sensed observations of the type described 
here in this paper it is not clear how this could be known. 

In Figure 3 we plot contours of ^/VL calculated from 
equation 7 for different stages of the shock wave evolu- 
tion in K — 9bn space, where k, is the index of the kappa 
distribution of the upstream ions, and Obn is the shock 
obhquity. AU panels assume T = 2 x 10^ K (173 eV), and 
a constant shock speed of 2100 km s~^. Panel 2a gives 
contours at 1.8 Rq, with Ma — 2 and va = 1050 km 
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s^^. Panel 2b gives contours at 2.1 Rq and Ma — 3, 
VA = 700 km s^^; 2c gives contours at 2.5 Rq and 
Ma ~ va = 525 km s~-^; and 2d gives contours at 
3.5 Rq and Ma — 8. The contour where the growth 
rate is zero is given by the thick dashed line. Other con- 
tours are given in units of the ion cyclotron frequency. 
In all cases, the requirements on k are least stringent 
at 9bn — 0°, and at the higher Ma shocks with higher 
plasma /3. At 1.8 Rq, we require k ~ 3.0 for paral- 
lel injection, rising to 4.6 at 2.1 Rq, 6.2 at 2.5 Rq, 
and 12 at 3.5 Rq, all in approximate agreement with 
equation 8. Early in their evolution, CME shocks are 
often expected to be quasi-perpendicular in geometry, 
and at 45° we require k ~ 2, 3.4, 4.5, and 9.5 respec- 
tively. In panels (c) and (d) , some of the contours show 
increasing k with increasing 6b„. We have verified that 
this behavior arises from our approximation concerning 
the value of fl/ui = max {Ma cos 6bm 5), and that the 
contour where the growth rate is zero is relatively un- 
affected by this. In reality, lu should be a spectrum of 
waves, and not a single value. The maximum value, 5, 
is motivated by the discussion of the Ma = 3.1 shock in 
|G_argate & Spitkovsky (2012), and should probably in- 
crease with increasing Ma, although it is kept constant 
here. 

The model of iGahnskv fc Shevchenkol (|20Tll ) assumes 
that waves generated ahead of the shock reach a steady 
state. This requires growth and damping rates of or- 
der n. This can be seen to be satisfied in Figure 3 as 
K — > 1.5 out to 9hn ~ 45° (t he obliquity considered by 
iGalinskv fc Shevchenko|[20Tll ) in panels c and d, but less 
so where Ma is smaller. This highlights the seed particle 
problem. Acceleration is inferred to commence at low 
coronal heights, but here shocks are generally of too low 
Ma to self inject particles from the upstream medium. A 
seed particle population is necessary at these altitudes. 

3.3. Remote Detection of Suprathermal Seed Particles 

We discuss the observation of such a suprathermal ion 
population directly by its effect on the profile of reso- 
nantly scattered Lyman a from the extended corona. 
Neutral hydrogen atoms remain coupled to the ion- 
ized component so long as the charge exchange rate 
^ niUcxV<Pv < u/r, where Ui is the proton density, acx, 
is the charge exchange cross section (e.g. Kadyrov et al. 
2006), V is the proton thermal velocity, while on the right 
hand side u/r expresses the solar wind expansion veloc- 
ity divided by the heliocentric radius. This translates to 
rii > 10^ cm~^, or a r < ARq. A population of seed par- 
ticles should reveal itself through extended wings on the 
Lyman a profile. Following arguments above, we model 
the distribution of coronal protons and neutral hydro- 
gen as a kappa distribution, and compute the profile of 
scattered disk Ly-a radiation to be observed with a suit- 
able instrument . We follow the formalism presented in 
ICranmeU (jl998[ ) , remembering the difference between his 
and our definitions of k. 

We model the Ly-a disk profile by 



Hz, and the central self-absorption. The ion density in 
the extended corona is given by 
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which matches Figure 2 in iCranmeii ()1998[ ). The density 
of neutral H is based on this equation, modified by the 
ratio of the radiative recombination rate (from iSeatonI 
1959fl and the electron collisional ionization rate (from 
Scholz fc Walter^ [1991 . The expansion velocity profile 
of the slow speed solar wind is taken to be 



20 + 300 logio (^/2^o) kms" 
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'l -0.95 exp (-10-23 (iy- i/o)^)) ergf 
which captures the Lorentzian wings with Ai/ = 4 x 10^^ 



which corresponds to Figure 8 in iCranmeri (pOOl ). 

In Figure 4, we show radiation transfer calculations 
following the formalism in Cranmer (1998), for the scat- 
tering of disk Lyman a radiation by hydrogen atoms in 
an isotropic "Kappa" distribution in the corona. The 
panels a, b, c, and d correspond to the profiles required 
for injection at Otn = 45° at 1.8, 2.2, 2.5, and 3.5 Rq, 
i.e. with the required values of k of 2, 3.4, 4.5, and 
9.5. In each case, the thick dashed histogram shows the 
kappa distribution profile, and the solid thin histogram 
shows the profile resulting from resonant scattering in 
a Maxwellian with temperature 2 x 10^ K (173 eV) for 
comparison. The thick solid histogram shows the profile 
resulting from having 0.1 of coronal H atoms in a kappa 
distribution and the remaining 0.9 in a Maxwellian. A 
Maxwellian distribution (where k — ^ oo) makes no con- 
tribution to growth or damping rates, see equation 8. 
The thin dotted curve shows the disk Ly a profile, to 
be read in ergs cm^^s^^sr^^A^^. The histograms are to 
be read in counts per bin, where the bins are 0.1 Awide. 
We have assumed a 10^ s integration time with a pixel 
size of 1 arcmin^, and an effective ar ea varying in the 
range 1-2 cm^ according to data in lOsantowski et ahl 
()1991t ). Any instrumental linewidth is neglected com- 
pared to the Doppler broadened coronal linewidth. The 
only noise included is that due to Poisson statistics in an 
assumed photon counting detector. Geocoronal absorp- 
tion is not include d, but would only affe ct the signal in 
one spectral pixel ([Meier fc Mangelll970D . Seed particles 
energized by prior shock passage might be expected to fill 
the observed volume to a greater extent than would seed 
particles associated with a reconnection current sheet or 
outflow shock. As discussed below this last possibility is 
perhaps the most plausible scenario for the generation of 
seed particles, and so one should concentrate on the solid 
histograms showing the 10:90 mix of seed particles and 
ambient plasma. Such seed particle distributions appear 
to be detectable out to 2.5 Rq, bearing in mind that 
longer integration times and larger spatial pixels are def- 
initely feasible. Since seed particles apparently have an 
effect up to 24 hrs after a prior CME (Gopalswamy et 
al. 2004), integration times an order of magnitude larger 
are definitely feasible, even with an observing duty cycle 
(lO^f 1/3 - 1/2. Thus not only would the detection of a 
suprathermal population be possible, but also its char- 
^A acterization in terms of density, k, and the likely ability 
to produce an SEP event when combined with a CME 
shock of appropriate speed. 
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4. ORIGINS OF CORONAL SUPRATHERMAL PROTONS 

So far we have discussed the existence of seed par- 
ticles from the point of view of the effects they would 
have on SEP production, and not said much about how 
they might be produced. A number of different ideas 
have been described in the literature. For the most part, 
these scenarios imply different expectations in the spec- 
tral, temporal, and spatial distributions of the suprather- 
mal protons and the Ly a spectra that they generate. We 
briefly outline here some of these scenarios and indicate 
how observations could distinguish among them. 

4.1. Suprathermal Protons from Preceding CMEs 

The result of Gopalswamy et al. (2004) suggests that 
previous solar activity (flare/CME) may prime the so- 
lar corona and wind with seed particles which can then 
be accelerated as SEPs by a second CME. This scenario 
should show significant differences in the coronal Ly a 
profile before and after an eruptive event and would per- 
haps be the easiest signature to distinguish should it be 
present. The seed particles may be the result of an ear- 
lier CME, but they are more likely SEPs associated with 
an earlier associated "impulsive" flare (Reames 1995). 
While the ambient solar wind composition is simply that 
of the corona, exhibiting the well known First Ionization 
Potential (FIP) enhancements of ions with low first ion- 
ization potential (e.g. Feldman & Laming 2000), impul- 
sive flare SEPs have a high Fe/0 abundance (Tylka & Lee 
2006a) suggesting that they are not ambient solar wind 
or coronal particles energized by a shock. More probably 
they are accelerated via reconnection at a flare-associated 
current sheet, thought to be a plausible site for the ori- 
gin of the abundance anomalies (Drake et al. 2009). The 
high Fc/0 particles also have high charge states which 
appear difficult to achieve simply by particle stripping 
during the acceleration process (Kocharov et al. 2000; 
Barghouty & Mewaldt 2000). Impulsive flare SEPs also 
exhibit increasing abundance enhancement over coronal 
values with increasing element mass (or possibly decreas- 
ing charge to mass). It is now known that this anomaly 
continues through the periodic table to the highest mass 
elements observed. Unlike the FIP effect, where an un- 
derstanding of the phenomenon is established and only 
details remain to be worked out (e.g. Laming 2004, 2009, 
2012), theories of impulsive flare abundances are much 
less secure. It appears that fractionation at reconnecting 
current sheets, either in the corona (Drake et al. 2009) 
or in the chromosphere (Arge & MuUan 1998), may yield 
the increasing enhancement with ion mass. 

Whether reconnection can supply particles of sufficient 
energies to be SEP seeds is another question. Generally, 
the reconnection exhaust is an outflow at the local Alfven 
speed, which will not meet perpendicular shock injection 
requirements. iGordovskvv et al] (|2010al lH) demonstrate 
using a test particle approach that some particles may be 
accelerated to much higher energies in the electric flelds 
of the reconnection current sheet. They find energetic 
protons in power law distributions with energy propor- 
tional to E^^^^^-^ , corresponding to k = 2 — 2.5. Alter- 
natively, magnetized particles with reconnection current 
sheets subject to the plasmoid instability may be further 
Fermi accelerated by the contraction of magnetic islands 
(Drake et al. 2006). Although this appears more plau- 



sible for electron acceleration than for ions (e.g. Drake 
et al. 2009), since to participate particle s need to be ini- 
tially super- Alfvenic. lDrake et all (|2012[ ) show from con- 
sidering a Fokker-Planck model of the process that an 
E~^'^ (corresponding to k = 2.5) spectrum is generally 
obtained, at least in the limit that the acceleration time 
is short compared to the loss time. 

Another plausible scenario appears to be ion acceler- 
ation at a fast mode shock occurring as the reconnec- 
tion outflow encounters denser plasma (e.g. Selkowitz 
& Blackman 2007; Mann et al. 2006), or at fast mode 
shocks produce d in the exhaust by an u nsteady recon- 
nection process ()Tanuma fc Shibatall2005D . Most authors 
consider electron acceleration by this means, with a view 
to modeling hard x-ray emissions observed during flares. 
We will discuss this in more detail below, but ion acceler- 
ation should also occur her e. Such a shock typica lly has 
a compression ratio r 2. ([Workman et al.l[2011[ ) and is 
quasi-perpendicular. iPark et al." f2012fl make the point 
that such a shock is typically high beta, because magnetic 
field has been destroyed in the upstream region by the 
reconnection process and converted to heat and kinetic 
energy. Thus the particle injection problem encountered 
for such a shock in low beta plasma (the subject of this 
paper for CME shocks) does not ap ply for the recon- 
nection outfiow shock. Figure 6f in (Park et all pOll 
shows the downstream ion distribution obeying a power 
law ^ E^"^ up to about 20 keV. This corresponds to 
and K c± 3, the values t o be e xpected for a shock with 
r = 2. iGuo fc Giacalond ()2012[ ) find even steeper proton 
spectra. With reference to Figure 2 here, the ion spectra 
associated with acceleration at a reconnection outflow 
shock appear to be too soft to provide the required seed 
particle distribution, particularly at the lower Ma val- 
ues found closer to the sun. Particles accelerated at an 
outflow shock are also not expected to be confined to the 
reconnection current sheet. This matches with the ob- 
served onset of ion acceleration at CME shocks occurring 
over a large area of the shock surface. 

The possibility also remains that seed particles for 
one CME could be produced at the forward shock of 
the preceding event. This would not easily produce the 
abundance variations with energy described above, un- 
less flare remnant particles from the reconnection site can 
be injected into the CME shock; unlikely since reconnec- 
tion exhausts move only at the ambient Alfven speed and 
ought not to be able to catch up with a shock. It is also 
not guaranteed to generate hard enough suprathermal 
particle spectra to excite upstream waves, similarly to 
the reconnection outflow shock discussed above. Another 
theoretical argument against producing seed particles at 
the shock itself is that shock acceleration times are in 
general short compared to the shock evolution time (e.g. 
Zank et al. 2006, Ng & Reames 2008). Given the dif- 
fusion coefficients, Hxx, of Zank et al. the acceleration 
timescale, Tacc oc t^xx I vly^^^f. is of order minutes to hours 
for a perpendicular shock. This is ample time to produce 
SEPs from a single shock. The fundamental obstacle to 
shock acceleration appears to be the injection of parti- 
cles. If a shock were easily able to do that, the survey 
of Gopalswamy et al. (2004) would have turned out very 
differently. 
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4.2. Suprathermal Protons from Ubiquitous 
Reconnection Processes 

In the modeling of Tylka & Lee (2006), small impul- 
sive SEP events, with their distinctive abundance distor- 
tions, were used to model the "flare" component of the 
suprathermal seed population. But impulsive events are 
not necessarily associated with a visible X-ray flare. In- 
stead, they are produced by reconnection processes in the 
corona, which may or may not be sufficiently energetic to 
accelerate electrons and generate X-rays. It is therefore 
possible that such reconnection activity goes on in the 
corona at nearly continuous level, generating suprather- 
mals but at levels that are far below detectability at in- 
terplanetary spacecraft. In this case, we would expect to 
see ubiquitous and omnipresent evidence for suprather- 
mal protons, perhaps mod ulated in intensity by the gen- 
eral level of flare activity. iMasson et al.l (|2013t ) consider 
the release of solar-flare accelerated particles onto open 
solar wind field lines as a fiux rope erupts and undergoes 
interchange reconnection with neighboring field lines. 

4.3. Suprathermal Protons from the Decay of 
Flare- Generated Neutrons 

Feldman et al. (2010) reported a surprisingly large flu- 
ence of solar neutrons observed by the Messenger space- 
craft at 0.48 AU following an M2 flare. Solar neutrons 
are produced when energetic particles accelerated in a 
flare loop travel downward and impact the chromosphere, 
with the resulting nuclear collisions liberating neutrons, 
some of which travel outward from the Sun. Feldman et 
al. also suggested that lower-energy neutrons that decay 
into protons while traversing the corona might provide a 
pool of suprathermal seed protons for shocks. Although 
there are concerns about background-contamination in 
the Messenger measurements (Share et al. 2011), the 
notion of seed protons from the decay of flare-generated 
neutrons may nevertheless be viable. If this process were 
to occur, we should expect the enhanced suprathermal 
proton population to appear within a few hours of the 
flare. The spatial distribution would probably be sym- 
metric about the normal to the flare location, with a 
neutron energy spectrum (and subsequent decay-proton 
energy spectra) dependent on the angle to the normal. If 
gamma-ray observations (from RHESSI or Fermi) were 
also available, modeling should be able to account for the 
observed radial distribution and energy spectrum of the 
suprathermal protons (Murphy et al. 2007, 2012). 

4.4. Magnetic Pumping Model 

One of the most intriguing discoveries in recent years 
has been the ubiquitous presence of suprathermal tails 
on the momentum distributions of ions in the solar wind, 
with a common spectral shape of in momentum. The 
common spectral shape occurs in the quiet solar wind far 
from shocks, in disturbed conditions downstream from 
shocks, and throughout the heliosheath. In the Fisk & 
Gloeckler model (2007, 2008; see also Fisk, Gloeckler & 
Schwadron 2010) this suprathermal tail is accelerated 
from the core distributions in compressive turbulence 
that is thermally isolated. It is further suggested that 
this mechanism may also serve to generate suprathermal 
ions in the corona, thereby providing seed particles for 
shock acceleration. If this model is correct, ubiquitous 



and omnipresent wings on the Ly a line profile should be 
observable, although the magnitude may vary with local 
turbulence levels. In addition, the p~^ spectrum will cor- 
respond to a distinctive spectral shape in the Ly a wings. 
It is also sufficiently hard to provide seed particles able 
to generate waves upstream of shocks, out to 6bn — 30° 
for Ma ^2 and Obn ^ 55° for = 3. 

4.5. Observing Strategy 

Gopalswamy et al. (2001) and Mann et al. (2003) 
model the Alfven speed above an active region by assum- 
ing the active region magnetic field to be a dipole super- 
imposed on the large scale solar magnetic field. Depend- 
ing on the CME speed, a shock is expected to form in the 
region 1.2-3 Rs- Above the first critical Alfven Mach 
number (Edmiston & Kennel 1984), a number of order 
1.5-2 depending on shock geometry, the shock must ac- 
celerate particles. Slit positions chosen below and above 
the height where many propagating disturbances associ- 
ated with CMEs are expected to steepen into a shock and 
begin to accelerate particles would allow discrimination 
between different production scenarios. Taking these slit 
positions at 1.5 and 3.0 Rq, if the CME shock produces 
seed particles we expect to see a signal first in the 3.0 
Rq slit. If the reconnection current sheet is responsi- 
ble, the 1.5 Rq slit should show the clearest signal. The 
calculations in Figure 4 are based on an assumed 10^ s 
integration in an arcmin^ pixel, which corresponds to the 
time taken for a 1000 km s~^ CME to travel 1.5 Rq. One 
could clearly use shorter integration times during a CME 
event to improve the discrimination of the timing of the 
appearance of Ly a wings. Figure 5 replots some of the 
examples from Figure 4 at shorter integration times, 10 
s and 100s. Signatures of seed particles are still visible 
at these shorter cadences. 

Such timing discrimination would also be important 
for inferring the presence of seed particles due to fiare- 
generated neutrons. The other two scenarios, suprather- 
mals from ubiquitous reconnection events and the p~^ 
proton spectrum resulting from thermally isolated com- 
pressive turbulence have no specific time dependence, 
but would have distinctive spectral signatures to observe. 
The observing strategy would probably involve long inte- 
gration times during quiet times to establish the possible 
existence of a seed particle population and to characterize 
its properties (i.e. value of kappa and number density). 
Given this information, one could predict the CME shock 
velocity required in order for an SEP hazard to develop, 
using the analysis presented above in section 3.1. During 
CME events observations should switch to shorter inte- 
gration times to take advantage of the relative timing of 
the appearance of suprathermal ions (if any) in the re- 
gions of the corona viewed by the different spectrometer 
slits. 

We end this section by mentioning another require- 
ment: seed particles produced by a prior CME must re- 
main with sufficient number density in the low corona to 
affect the evolution of a second CME for a period of up 
to about 24 hours, a timescale about an order of mag- 
nitude longer than that expected if seed particles were 
swept out with the solar wind. Here we estimate the dif- 
fusion coefficient required to restrict the outward motion 
of energetic particles. 

The relevant diffusion coefficient is (e.g; 



8 

IBlandford fc EichledHOSTl: IRakowski et alJl2008[) 



where = qB/mc is the particle cyclotron frequency, 
U {k\\ = r2/w||) = p6v\^J1 is the wave energy density 

at parallel wavenumber fcy = r^/wy, and other symbols 
have their usual meanings. Taking B = IG, a num- 
ber density of n = 10^ protons cm^'^, and ^ w/2, 
I? ~ 3 X V^^'^v'^ 1 8v\^^ cm^s"-^. Putting the diffusion 

timescale t = jD with lengthscale I ~ li?©, we find 
t ~ 2 X lO^Jw^ii /w^ Setting t ~ 8 x 10"* s, 8vk^^ /v ~ 

2 X 10~^. So trapping a 10 keV proton {v ^ 10^ cm s~^) 
requires 5vk^^ — 20 km s~^. This is of similar order of 
magnitude, but smaller than the thermal line broadening 
at coronal temperatures, and so probably not unambigu- 
ously detectable as a signal of seed particles. 

5. SUPRATHERMAL ELECTRONS 

Many of the mechanisms discussed above also acceler- 
ate electrons. While the necessity of electron seeds for 
further shock acceleration is controversial, the observa- 
tion and characterization of any suprathermal electron 
distribution may further aid in understanding the pro- 
cesses responsible for any suprathermal ion population. 

Recent observations (e.g. Krucker et al. 2010) have 
indicated that the energy in accelerated electrons dur- 
ing a flare can be similar to that in the magnetic field, 
and that the whole population of electrons is energized. 
This argues strongly in favor of electron acceleration by 
a Fermi process in contracting magnetic islands in a re- 
connecting current sheet(Drake et al. 2006, where such 
a degree of electron energization was actually predicted) 
rather than acceleration at an outflow shock (where only 
a small fraction of the particles would be expected to 
be accelerated). In the limit of strong magnetic field, 
or electron plasma /3 — 0, a v"^ electron distribution 
is predicted, becoming a steeper power law when the 
back-reaction of the electron pressure on the acceleration 
process is included. Consequently reconnection in lower 
electron plasma /? conditions is expected to put more of 
its released magnetic energy into accelerated electrons. 

In Figure 6 we show various illustrations of how heated 
electrons distort the extended wings of the Lyman a pro- 
file through Thomson scattering. Panels a and b show 
the profiles of a Maxwellian at temperature 2 x 10^ K 
(173 eV), and kappa distributions with k — 2.5 and 
5, intended to match the r ange of electron spe ctra pro- 
duced by the simulat i ons in iDrake et al.l ()2006[ ) , and ob- 
served by lOka et all ()2012[ ) . Panel a gives the case at 
1.8 i?0, and panel b at 2.2 Rq. It can be seen that 
K = 5 in Thomson scattering is indistinguishable from 
a Maxwellian. The k = 2.4 profile however is distinct, 
and potentially detectable. We have assumed a 10^ s in- 
tegration time, with emission concentrated in a 0.1 xl 
arcmin^ pixel, assuming that we would be observing elec- 
trons heated in a reconnection current sheet of about 0.1 
arcmin width. The histograms have a 1 A bin size. Pan- 
els c and d show similar calculations for the same loca- 
tions, but with Maxwellian electrons with temperatures 
of 2 X 10^, 4 X 10^, and 8 x 10^ K. These different pro- 



files are more distinguishable from each other, and offer 
further diagnostics of reconnection processes. 

6. SUMMARY 

SEPs are a major radiation hazard for both robotic 
and manned missions. The episodic nature of the SEP 
hazard presents engineering challenges different from 
those of the other ionizing-radiation populations found 
in space. Discoveries about the nature, variability, and 
origin of coronal suprathermal protons would be signifi- 
cant breakthroughs in understanding the provenance of 
SEPs. We have argued based on prior observational ev- 
idence and theoretical considerations that the presence 
of a suprathermal seed particle population is an impor- 
tant factor, and maybe the crucial factor, in determining 
whether a CME becomes "SEP effective" . The observa- 
tion of such a population, through for example the line 
profile of Lyman a emitted from the extended corona, 
would be a major step forward in our understanding. 
We further anticipate that similar instrumentation, ulti- 
mately deployed as a monitoring mission, would greatly 
enhance operational forecasts of SEP events, and also 
provide notice of "all-clear" periods, through the pres- 
ence or absence of seed particles as revealed by spec- 
troscopy of the H I Lyman a line profile. 

Given that Cycle 24 is presently predicted to have a 
sunspot maximum in June- July 2013 that will be smaller 
than any other observed during the space age, it is diffi- 
cult to extrapolate from previous experience to say how 
many SEP events would be observed in the next few 
years. However, it is worth noting that some of the 
largest SEP events in the historical record, including Au- 
gust 1972, October 2003, January 2005, and December 
2006, occurred in the declining phase of the cycle. An 
experiment launched in the next few years would be ex- 
pected to be on-orbit while the new Solar Orbiter and So- 
lar Probe missions are operating. These missions, com- 
bined with others likely to be operational in an extended- 
mission mode (SDO, Hinode) will provide data on the so- 
lar and inner-heliospheric context of such observations. 

Some of our work al so has astrophysical applications. 
iRavmond et al.l ()2010t ) observe a kappa distribution in 
emission in H a from the forward shock wave in Tycho's 
supernova remnant. This arises as neutral H in the up- 
stream medium passed through the shock, and is excited 
to emit H a before it is ionized postshoc k, and thus re- 
tains its preshock distribution function. iRavmond et al.l 
use a different definition of k, so their quoted value 
of K ~ 2 corresponds to k ~ 3 in our work. At this value, 
we estimate a fraction of order 10^'^ — 10^^ of the shocked 
plasma particles are able to join a seed particle popula- 
tion. Raymond et al. (2010) express some concern that 
such a value of n implies very efficient particle accelera- 
tion. Considerations of seed particle injection outlined in 
this paper suggest that this need not be the case, espe- 
cially at oblique shocks close to the perpendicular limit, 
though of course the Alfven Mach number is likely much 
higher. 
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Fig. 1. — Evolution over the first ten minutes of the proton spectrum at the shock (based on Ng &; Reames 2008) as a CME moves 
outward from ~ 3Rq at 2000 km s~^ for two different levels of input suprathermal seed protons (first panel). Different colors in the second 
and third panels correspond to different elapsed times in seconds since the start of particle acceleration, as given in the labels. These 
calculations take into account growth of proton-amplified Alfven waves, which is a non-linear process. Accordingly, the simulation results 
are in real units and are not adjusted by any arbitrary normalization factors. 




Fig. 2. — The "kappa" distribution for various values of the k parameter. Curves are shown for k = 2.5, 5, 10, 20, and also the limiting 
case K — > oo where a Maxwellian is obtained. The velocity axis is given in units of the thermal speed vxh entering equation 1. 
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Fig. 3. — Contours of upstream wave growth rate given by equation 7 in k — 0i,„ space, where k is the index of the kappa distribution 
of the upstream ions, and Sf,^ is the shock obhquity. All panels assume T = 2 X 10 K (173 eV), and give in order the cases Ma = 2, 
Va = 1050 km s^i; Ma = 3, Va = 700 km s^^; Ma = 4, Va = 525 km s^^; Ma = 8, Va = 262.5 km s"!; corresponding to parameters for 
the shock of the 2005 January 20 event at 1.8, 2.1, 2.5 and 3.5 Rq respectively. The contour where the growth rate is zero is given by the 
thick dashed line. Other contours are given in units of the ion cyclotron frequency. In all cases, the requirements on k are least stringent 
^■t ^6n = 0°i ^■iid at the higher Ma shocks with higher plasma /3. At 1.8 Rq, we require k ~ 3.0 for parallel injection, rising to 4.6 at 2.1 
Rq, 6.2 at 2.5 Rq, and 12 at 3.5 Rq. At these altitudes, injection can extend out to about 35°, 60°, 65°, and 80° respectively. 
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Fk;. 4. — Simulated Lyman-o profiles corresponding to the K-distributions required for injection into the parallel shock regions in the 
panels of Figure 3. The radiant profile is indicated as a dotted line, to be read on the intensity axis in units ergs cm^^s^^sr~^A~^ . The 
dashed histogram shows the result of scattering of this radiation in a K-distribution of hydrogen atoms with value 2, 3.4, 4.5, and 9.5 
respectively. The thin histogram shows the result in each case for a Maxwellian, with the thicker histogram showing an intermediate result 
with 90% of the hydrogen atoms along the line of sight in a Maxwellian distribution, and 10% in the K-distribution. We have assumed 
Poisson statistics in a photon counting detector, with eflfective area of ~ 1 cm^, and 10'^ s integration time. All histograms give the counts 
in 0.1 A bins, in a spatial pixel of 1 arcmin^. Such seed particle distributions are easily detectable out to 2.5 Rq, bearing in mind that 
longer integration times and larger spatial pixels are definitely feasible. 
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Fig. 5. — Simulated Lyman-a profiles from Figure 4a integrated for 10 and 100 s (panels a and b), and from Figure 4b and 4c integrated 
for 100 s (panels c and d). The extended wings are still detectable at these shorter integration times. 
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Fig. 6. — Simulated extended Lyman-a profiles showing the effect of Thomson scattering by hot electrons. Panels a and b show the 
results at 1.8 Rq and 2.2 Rq at T = 2 x 10 K (173 eV), for a MaxwcUian (thin solid histogram), k = 5 (thick solid histogram) and k = 2.4 
(thick dashed histogram). Panels c and d show results at 1.8 Rq and 2.2 Rq at T = 2 x 10^ K (thin solid histogram), T = 4 x 10^ K (350 
eV; thick solid histogram), and T = 8 x 10® K (700 eV; thick dashed histogram). The assumed integration times is 10^ s, in a pixel of size 
0.1 x 1 arcmin^, determined by the angular width of a reconnection current sheet. The histograms have 1 A bins. 



